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A mechanism of self-heating of solar corona is pointed out. It is shown that the free 
energy available in the form of sheared flows gives rise to unstable electrostatic waves 
which accelerate the particles and heat them. The electrostatic perturbations take 
place through two processes (a) by purely growing sheared flow-driven instability and 
(b) by sheared flow-driven drift waves. These processes occur throughout the corona 
and hence the self-hcating is very important in this plasma. These instabilities 
can give rise to local electrostatic potentials <^ of the order of about 100 volts or 
less within 3 x 10~^ to a few seconds time if the initial perturbation is assumed to 
be about one percent that is ^ ~ 10~^. The components of wave lengths in the 
direction perpendicular to external magnetic field Bq vary from about 10m to Im. 
The purely growing instability creates electrostatic fields by sheared flows even if 
the density gradient does not exist whereas the density gradient is crucial for the 
concurrence of drift wave instability. 

Subject headings: Sun: self-heating of corona, sheared flow-driven instability, drift 
waves. 

PACS numbers: 

I. INTRODUCTION 



Several theoretical models have been presented to explain the cause of solar 
coronal heating (Aschwanden 2001; Mandrini, Demoulin, h Klimchuk 2000). 
It is a well-known counter institutive fact that the corona is 200 times hotter 
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than the chromosphere while it is a rarefied collision-less plasma with electron 
temperature Tg ~ 10^ (K). There are many puzzles in the problem of coronal 
heating. The temperature rises by two orders of magnitude from upper chro- 
mosphere to lower corona through a transition layer of only about 500 km 
(Priest 1982; Narain; Ulmsch- neider &; Klimchuk 2006). The proposed wave 
heating mechanisms assume that the waves originate in the lower regions and 
deposit their energy in the corona. Observational data of Skylab, Yohkoh, 
SOHO and TRACE has indicated that the entire corona is filled with open 
and closed magnetic field lines and only a subset of it is loaded with hot 
plasma at a given time (Litwin & Rosner 1993; Hara et al. 1992; Moses et 
al. 1997; Schrijver et al. 1999). The coronal loops have also been considered 
to be responsible for localized heating. 

The bright coronal loops have higher densities compared to the ambient faint 
coronal plasma which indicates that the heated plasma originates from the 
dense chromosphere. The coronal heating by acoustic waves originating from 
global oscillations have been ruled out (Aschwanden 2001). Alfven waves 
have been considered to be the best candidate for carrying adequate energy 
fluxes from chromosphere to corona (Hollweg & Sterling 1994). A great deal 
of work on Alfven wave heating of the corona has appeared in literature (Ion- 
son 1983; Mok 1987; Steinolfson & Davila 1993; Ofman, Davila & Steinolfson 
1995; Helberstadt k Goedbloed 1995; Ruderman et al. 1997). 
Any unstable wave in this plasma can cause coronal heating through damping 
either by wave-particle interaction or by other linear and nonlinear mecha- 
nisms. Most of the theories pertaining to wave heating in corona are based 
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on magnetohydrodynamics (MHD) and hence electrostatic drift-waves have 
not been investigated while in fact the system itself is highly inhomogeneous. 
There are many suggested mechanisms for heating in which the main energy 
source is provided to corona by other regions. 

We point out an important energy source available within the corona to ex- 
cite short scale electrostatic perturbations and that are the sheared flows of 
plasma streams and rivers. These electrostatic perturbations can be associ- 
ated with either the purely unstable mode (D'Angelo 1965) or the drift- waves 
of two fluid plasma depending upon the scales of density inhomogeneity and 
the wavelengths. Vranjes and Poedts (2009a) have proposed a new paradigm 
of solar coronal heating by drift waves. They have used the results of ki- 
netic theory which predicts that drift wave is a universally unstable mode 
and density gradient is the source of its instability which is also the cause of 
its existence. It has been proposed that the drift waves can heat the corona 
through two possible ways, one due to the Landau damping effects in the 
direction parallel to the magnetic field, and another one, the stochastic heat- 
ing in the perpendicular direction. It has been confirmed (De Pontieu et al. 
2007) that the solar atmosphere is highly structured and has inhomogeneous 
density filaments of various sizes. Therefore drift waves in corona may have 
several different scales of wavelengths. 

Vranjes and Poedts (2009a) assume external magnetic field is along z-axis 
while the density gradient is along x-axis and wave propagates under local 
approximation in yz-plane with perpendicular wavelength Xy = ^ P^^"" 
allel wavelength = ^ while kz « ky (where ky and kz are respective 
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wave numbers). The drift waves have been considered with \y ~ 0.5m while 

-1 

km) is the density gradient 
scale length. 

We focus our attention on the fact that sheared flows are omnipresent in the 
corona. Thus a large source of energy for the excitation of short scale elec- 
trostatic perturbations exist within the corona. Two types of electrostatic 
instabilities are expected to occur throughout the corona continually: 

I. Purely growing sheared flow-driven instability in doppler shifted frame 
(D'Angelo 1965) which exist even if the plasma density is uniform 

II. Drift wave instability which needs plasma density gradient for it's existence 
(Kadomtsev h Timofeev (1963)) and sheared flows for instability (Saleem, 
Vranjes k Poedts (2007)) 

We investigate the drift waves having relatively longer wavelengths com- 
pared to ion Larmor radius so that the fluid model can be justified. More- 
over, the electrons are assumed to follow Boltsmann distribution. Therefore 
frequencies Ur of drift waves must fulfill the condition Ur « J^tekz where 
= (Te/rrie)^/^ is the electron thermal velocity and kz is the wave num- 
ber parallel to external magnetic field. The drift waves and purely growing 
D'Angelo mode have shorter wavelengths compared to Alfven waves of MHD 
model, in general. The electrostatic perturbations driven by sheared flows 
can transfer their energy to plasma particles. For the case of drift waves 
a detailed picture has been presented by Vranjes & Peodts (2009a). The 
drift dissipative instabilities have also been discussed by Vranjes & Peodts 
(2009b). 



= 20A;m and L„ = 100m (where L„ = 



no dx 
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II. THEORETICAL MODEL 

Let us consider the collisionless coronal plasma consisting of two types of 
ions j = a, b, where 'a' represents Hydrogen ions and 'b' represents Helium 
ions. We choose the ambient magnetic field as Bq = Bqz =constant, Vn^o = 

dujo 



— X , and the shear flow as i'o{x) = Vjo{x)z = Vq{x)z where vq is the same 
for all species; electrons and both type of ions. Each species has a zero-order 
diamagnetic drift vdjo = — (^^) V/nn^o x z where Tj is temperature and 
Qj is charge of the jth species of ions. For electrostatic perturbation, the 
momentum equation yields the perpendicular component of the velocity of 
jth species of ions as, 

v., = J^(E, X z) - ^(V/nu, X 2) - ^ (^) - j5:(a. + v,.V)v, x z 

= ve + yoj + vnj + Vpj (1) 

where ^^E,^Dj,^uj and v^j are the electric, diamagnetic, stress tensor and 
polarization drifts, respectively. The parallel component of momentum equa- 
tion gives, 

o ■ T- 

{dt + VjOzdz)vjzi + Vjzidxiyjozix) = —Ezi —dzUji (2) 

TTij TRjTljQ 

Here dx = and the subscripts naught (0) and one (1) denote zero order 
and linear quantities, respectively. The continuity equation for jth ions is, 

T 

^T^i^t + i^jOzdz)^'inji + rijodzVjzi = (3) 
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where Qj = and Vz^oz(a^) = +x 
the following relation 



dx 



In obtaining (3), we have used 



vo{x)dz 







QjBo rijo J 



The continuity equations for j = a, b can also be expressed as, 



(4) 



(5) 
(6) 



where 



^6 = (1 + PiTK)^t - ^tek' 



C.2 



+ Pbs^l^t + c-lgkyk 



2 , 2o2 ,2 7 7 1 ^^0 



2 l2 



Here we have defined v'^- = 

1 dn 



rife dx 



= vok 



rijo dx 



, and $ = ^. 



Poisson equation in this case can be written as, 



V.Ei = —{Ual + Ubl - Uel) 

Co 

and electrons (e) are assumed to follow the Boltzmann relation. 



(7) 



(8) 
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In steady state the relation ngo = nao + n^o holds. Equations (5 - 8) yields a 
fourth order dispersion relation as follows, 



L40i + LaQj + 12^1 + L^nl + LiO, + Lq = (9) 



where 



T A 2 2 2 I ^aO / 2 2 7 2\ I / 2 2 7 2\ 



r j ^aO 2 * , ^60 2 * 

L3 = - <^ — «6^a + — «a^6 
L ^eO ^eO J 

L2 = -A [a^iybrk, + (^h^aT^z) + — {a^QaCsakz 

_ 2 7.2 2 7.2 X , !^/^2 2 7.2 _ 2 7.2 2 7.2 x 

7- ^aO / 2 7 2 *\ , / 2 7 2 *\ 

neo rieo 

T _ a 2 2 7.2 2 7.2 _ ^ 2 7.2 2 7.2 _ ^„ ,.2 7.2 2 7.2 
^0 — ^aT'^z'^bT'^z ^ ya'^bT^z^as'^z ^ 9b'^aT'^zHs'^z 

neo neo 
Here g, = [^A, -l),Aj = i^,a| = (1 + p,VA^),AL = ^, and = 

(1 + ALfc')- 



III. APPLICATION TO CORONA 

Now we show that the theoretical model presented above is perfectly appli- 
cable to coronal plasma. Let us choose the parameters of the corona (Priest 
1982) as Ueo = 1015m~^, n^o = 0.9neo, and Tg = lO^K. Since the condition 
Te < Th < Tffe holds in this plasma (Hansteen, Leer, and Holtzer 1997), 
therefore we assume Ta = 2.5Te and = STg. 

First we show that the sheared flows of electron proton plasma of corona 
gives rise to purely growing and oscillatory drift wave instabilities. Then it 
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will be shown that the presence of second ion (Helium 10%) in this plasma 
modifies the growth rate and the real frequencies of the drift wave in differ- 
ent parameter regimes. If this 10% concentration of the ions in the plasma 
is neglected, then we may use ribQ = and T^o = in the equation (9). In 
the limit A|,g/c^ << 1, we have = 1, = 1, ^ = 1, u^j^ = and it gives 
^4 = plskl + L3 = -t^*, L2 = -nulr^kl + Qaclskl and Li = = Lq. In 
this case equation (9) reduces to 



(1 + pirK + pls^D^t - ^l^u: + Aak,ci, - [cl + vItW, = (10) 



This quadratic equation has two roots 

1 



2Af 



< ± <! + ^^ociM ( (1 + aa) - A, 



ky 



(^^a;)l,2 = 



1 



2Ao 

,2 l2 



Kf + 4Aoclkl ({1 + aa) - A^^ 



(11) 



where Aq = (1 + plr^y + Pas^y)- Following conditions should satisfy simulta- 



neously for the instability. 



{l + <Ja)]^<Aa 

Ky 



and 



Aa — 

kp- 



(12) 



(13) 



We shall assume that Knj has the same value for all species and that can be 
denoted by Kn- In Fig. 1, the growth rates of the shear fiow-driven elec- 
trostatic instability are plotted using equation (10) against different gradient 
scale lengths of the plasma flow parallel to external magnetic fleld in case of 
homogeneous density. The real frequency in laboratory frame is ujq. But in 
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the moving frame of plasma, it is a purely growing instability and therefore Ui 
has been plotted vs perpendicular wavenumber ky. We notice that the imag- 
inary frequency uji decreases corresponding to the same ky when the shear 
flow gradient scale length Ly = ^ increases. Thus the steeper gradients give 
rise to larger growth rate of the instability. 

In the presence of density gradient, the electrostatic drift wave can also be 
excited. Fig. 2 shows that the drift wave becomes unstable for 2 < /cy < 2.5 
for a very short range of wavelengths. Somewhere in between ky = 2 and 
ky = 2.5, the real frequency of drift wave uj* is larger than the growth rate uJi 
that is UJi < cj* and hence linear analysis is valid. The drift wave seems to be 
driven by the sheared-flow. The pattern of growth rate ui for 2 < ky < 2.5 
becomes different in Fig. 2 which has Kn 7^ ^ compared to Fig. 2 plotted for 
Kn = 0. For 2.5 < ky, the purely growing instability dominates because here 
Ur « UJi and hence it shows the shear flow-driven instability with local real 
frequency ujq = v^kz- Both uJi and ujr becomes larger for drift wave when the 
shear flow gradient is steeper as shown in Fig. 3 while the density gradient is 
kept constant = 1-9 x 10~^m~^. We notice that the growth rates and real 
frequency of the drift wave decrease corresponding to smaller ky and kz com- 
pared to Fig. 2 which is quite natural. In this case the instabilities require 
larger value of flow because c^g^l - term becomes smaller in condition (13). 
In Fig. 4, the efi"ect of the ration of ^ on instabilities is shown. In Fig. 5 
the effect of the presence of helium ions on real and imaginary frequencies is 
shown on the shear flow-driven instability (for k„ = 0) and on the drift wave 
instability for [tin 7^ 0). 
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IV. ROLE OF DISSIPATION 

The Solar corona is commonly assumed to be a collision-less plasma be- 
cause most of the wave studies deal with the relatively higher frequency 
Alfven waves using MHD equations or ion acoustic waves with the wave num- 
bers larger than our regime of parameters and hence the frequency becomes 
larger than the electron-ion collision frequency z/ei- The drift waves have been 
investigated in corona first time (to the best of authors knowledge) by Vran- 
jes and Poedts (2009a) and (2009b). The drift waves investigated by Vranjes 
& Poedts (2009a) have frequencies a;^- ~ a;^ ~ lO'^rad/s and i/gi — SOrad/s 
(Vranjes & Poedts 2009b). Therefore the ideal plasma approximation is valid. 
In the higher frequency {oj*a) regime it does not seem preferable to assume 
electrons to be inertia-less because the condition cj* << Vtekz may not remain 
valid. We choose the wave parameters such that the condition a;*, o;^ << Vtekz 
remains valid and electrons follow the Boltzmann distribution. But here we 
see another small effect namely the dissipation. Let us look at the Fig. 2(a) 
to analyze the role of drift wave in corona. This wave is stable for ky < 2 
and then for 2 < fcy < 3, the wave develops oji and for S < ky we can see that 
ujr < uJi and the shear flow-driven instability dominates. To understand the 
wave behaviour, we choose a value of ky in between 2 and 3. Let ky = 2.2m~^, 
then for k„ = 1.9 x 10~^m~^ we find a;* = DgKnky ~ 36rad/s which is equal 
to Vei- Therefore, dissipation can play some role in this frequency regime. 
It is important to note that the well-known drift dissipative instability (DDI) 
(Weiland 2000) does not become important even in this range of frequen- 
cies because a;* ^ Uei. The parallel momentum equation for electrons yields 
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(Weiland 2000) 

in the limit a;* << Uei- In ion continuity equation the hmit c^/c^ << a;* is 
used along with Ue — rii, to obtain linear dispersion relation for DI with the 
real frequency Urd and imaginary frequency ujid as (Weiland 2000), 



and 



^id - I -2T9 ^eiplskl (16) 



^2 y •^eit-as"'y 



In our case, the drift wave frequency can be of the order of Vei- Since we have 
shown that the drift wave becomes unstable due to sheared flow in the absence 
of electron-ion collisions, therefore in our parameter regime, the dissipation 
can just add its small effect to the already unstable perturbation. But the 
DDI is not applicable because the real frequency a;* is not much smaller than 



The Fig. 2(a) indicates that for small /Cy— values, the drift wave is stable 
and shear flow does not have an effect on it. So for longer wavelengths 
ky « 2, the drift dissipative instability in corona can take place. As an 
example, let us choose ky = O.lm^^ and take rest of the parameters to be 
the same i.e. Kn = 1.9 x 10~^, kz = lO'^^ky, Bq = 10~^T etc. Then we find 
cj* ~ 1.63, i^tekz — 36 while i^ei — 36. Then cj* ~ v^i « holds along 
with a;* << Therefore, the drift dissipative instability in corona gives 



rise to drift waves having very low frequency u* ~ IHz and relatively longer 
y k,. 



wavelength Ay ~ |^ ~ 60m 
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V. DISCUSSION 

It has been proposed that the large perturbed electrostatic fields are gen- 
erated throughout the solar corona, due to localized sheared plasma flows, 
which accelerate the particles and heat them. The corona is not a static ball 
of plasma, rather it has flows and gradients. Therefore, it has already been 
proposed that the free energy available in the form of density gradients can 
produce electrostatic drift waves (Vranjes & Poedts 2009a, 2009b). Using 
the results of kinetic theory, these authors have shown that the universally 
unstable drift waves can heat the coronal ions very efficiently. The density 
gradient in the direction perpendicular to the external magnetic field is the 
cause of the existence of these waves as well as it is the source for their insta- 
bility. The waves transfer their energy to plasma particles through Landau 
damping (the wave-particle interaction). This process cannot be studied us- 
ing fiuid models. 

But the two fiuid theory has predicted a Kelvin- Helmholtz type instability 
which takes place in plasmas because of sheared fiows (D'Angelo 1965). It 
shows that if both electrons and ions fiow with the same velocity along the 
external magnetic field Bq and there exists a gradient in fiow in the direction 
perpendicular to Bq, then the perturbed electrostatic field becomes unstable. 
These are purely growing fields in the frame of reference of the fiow wide spec- 
trum of wavelengths. But in laboratory frame these unstable perturbations 
have associated local real frequencies as ujq = vo{x)kz. The sheared fiows in 
corona can give rise to two types of electrostatic instabilities continuously 
throughout the corona. 
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1. Sheared flow-driven instability (D'Angelo 1965) even if the plasma den- 
sity is uniform {tin = 0) 

2. The drift wave, which exists if the density is non-uniform Kn 0, insta- 
bility due to sheared flow 

These instabilities have been investigated neglecting the eflFects of electron-ion 
collisions. For ky ~ 10m~^, we find that the sheared fiow-driven instability 
can create potential (/? ~ 68 volt in about growth-time Tg ~ 0.03s if at t = 
we assume ^ ~ 10~^. Note that Vranjes & Poedts (2009a) have estimated 
that large frequency drift waves {ujr ~ 2.5 x 10^) can give rise to this value of 
potential in 0.02 s. The components of wavelength chosen in the perpendic- 
ular and parallel directions are, respectively. Ay = 0.5m~^ and = 20/cm. 
Then they have vtekz = 1.2 x 10^. 

We have also shown that the drift waves having ujr < 50 rad/s satisfy 
ujr « ^tekz can become unstable due to sheared fiow. The collisional ef- 
fects have been neglected which can only modify the growth rates by small 
amounts. It is well-known that the electron-ion collisions can drive drift 
dissipative instability and it does not require plasma flow to occur. The 
collisions produce electrostatic drift waves having very small real frequency 
ujrd — 1.6 rad/sec corresponding to ky ~ 0.1m~^ and longer wavelengths 
Xy = Y — 60m. If the initial perturbation is assumed to be $0 = = 10" , 
then ip will take time ~ 15 minutes to grow up to = 86 volt. But smaller 
values of ip will be produced in much small times than 15 minutes. Thus we 
conclude that drift waves of different frequencies and wavelengths are pro- 
duced in the solar corona due to sheared flows and electron-ion collisions. 
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Fig . 1 

FIG. 1: (Color online) The growth rate (wj) is plotted against the perpendicular component of the 
wavenumber {ky) for different inverse velocity scalelengths; = ky /60 (solid curve), = fcy/lOO 
(dashed curve), and Ky = ky/200 (dotted curve) with rieo ~ nao ~ 10^^ m~^, Te = 10^ K, Ta = 2.5 
Tg, -Bo ^ 10^^ Tesla, vq = 10 km/s, Knj = 0, nf,o = 0, and k^ = 10^'^ky m^^. 

Thus electrostatic fields are almost omnipresent in the coronal plasma and 
continuous self-heating is taking place due to sheared flows and density gra- 
dients. The present investigation shows that even if the plasma density is 
uniform in a region, the electrostatic fields will be produced because corona 
is not static and sheared flows occur everywhere. Therefore, the sheared flow 
driven instability plays an important role in self-heating of the corona. 
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FIG. 2: (Color online) The real and imaginary frequencies \bJr (dashed curve), u;, (solid curve)] are 
plotted against the perpendicular component of the wavenumber {ky) for varying inverse velocity 
scalelength (a) = ky/60, and (b) = ky/100, taking vq = 10 km/s, k„ = 1.9 x 10^^ m~^, 
iT'bo = 0, and = lO^^/c^m^^. All other parameters are the same as in Fig.l. 



Captions 



Figure 6: (Color online) The imaginary and real frequencies (a) Ui and (b) 
LOr^ are plotted against the perpendicular component of the wavenumber {ky), 
respectively, for rifjQ = (dashed curve) and ribo = O.lngo (solid curve), with 
vq = 10 km/s, TT-eo ~ 10^^ m~^, Uao ~ 0.9neo and Bq ~ 10^^ Tesla, taking 
Kn = 0, Ky = and kz = 10~^/cym^^. 
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FIG. 3: (Color online) The real and imaginary frequencies [ujr (dashed curve), Ui (solid curve)] are 
plotted against the perpendicular component of the wavenumber (ky) for changing inverse velocity 
scalelength (a) k„ = ky/QO and (b) Ky = ky/100, keeping vq = 70 km/s, k„ = 1.9 x 10^'^ m~^, 
iT'bo = Oi Slid kz = lO^^kyin^^ . All other parameters are the same as in Fig.l. 
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FIG. 4: (Color online) The growth rate is plotted against the perpendicular component of 
the wavenumber (ky) and the inverse density scalelength respectively, varying the parallel 

component of the wavenumver (a) kz = 10~^A;j^m~^ (solid curve), k^ = 10~^A;j^m~"'^ (dashed curve) 
with fixed values of vq = 10 km/s, k„ = 1.9 x 10^^ m^-*^, and = ky/60, and the streaming velocity 
(b) Vq = 50 km/s (solid curve), vq = 70 km/s (dashed curve), with = ky/60, ky = 0.5 m"-'^, and 
kz = 10~^A;ym^^. All other parameters are the same as in Fig.l. 
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FIG. 5: (Color online) The imaginary and real frequencies (a) uji and (b) are plotted against 
the perpendicular component of the wavenumber (fcy), respectively, for n{,o = (dashed curve) and 
Tifoo = O.lneo (solid curve), with vq = 10 km/s, ngo ~ 10"*^^ m""^, ~ 0.9neo and Bq ~ 10~^ Tesla, 
taking k;„ = 1.9 x 10^'^m^-'^, = ky/60, and kz = W~'^kym~^ . 



Ofman, L., Davila, J.M. & Steinolfson, R.S. 1995, ApJ, 444, 471. 

Priest, E.R., 1982, Solar Magnetohydrodynamics (England, D. Reidel Publishing Company), 7. 
Ruderman, M.S., Berghmans, D. Goossens, M. k Poedts, S., 1997, A & A, 320, 305. 
Saleem, H., Vranjes, J., & Poedts, S., 2007, A&A, 471, 289. 
Schrijver, C.J., et al. 1999, Sol. Phys., 187, 261. 
Steinolfson, R.S. & Davila, J. 1993, ApJ, 415, 354. 
Vranjes, J. & Poedts, S. 2009a, EPL, 86, 39001. 



19 



Vranjes, J. & Poedts, S., 2009b, Mon. Not. R. Astron, Sol., 398, 918. 



